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NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL MEM0RANDUI4 1282 


THEORY OF FLAME PROPAGATION* 
By Y. B. Zeldovich 


The mechanism of flame propagation has been qualitatively 
formulated. In accordance with this formulation, the chemical 
reaction initiated in some layer brings about an increase in the 
temperature^ because of the heat conduction, the temperature is raised 
in the neighboring layer where in turn the chemical reaction is ini- 
tiated. In this manner the flame is propagated. 

Michelson (reference l) is credited with presenting the first 
reliable method for measuring the flame speed and for computing the 
temperature distribution in the preheated zone of burning gas. In 
various papers on the theory of the speed of combustion (Jouguet 
(reference 2), Nusselt (reference 3), Daniell (reference 4), and 
others), a solution is found only by making simplifying assumptions 
that do not correspond to the actual conditions. The comparatively 
unl-inown work of Taffanel (reference 5) is the most similar to present- 
day views. 

In the work of Jost and von Muffling (references 6 and 7), Sachsse 
(reference 8), and particularly Lewis and von Elbe (reference 9), a 
method is formulated that is satisfied by theory, that is', the dif- 
fusion as well as the heat conductivity is considered, the concept of 
ignition temperature is not included, and the continious dependence of 
the speed of reaction on the temperature and concentration is taken 
into account . 

This method was first carried out in the work of D. A. Frank- 
Ka.menetskii and the author (reference 10) , in which an approximate 
method is given for the computation of the flame speed xmder definite 
rational assumptions concerning the chemical reaction and the constants 
of the explosive mixture. 

The experimental work carried out at the Institute of Chemical 
Physics (reference ll) confirmed the validity of the theory. 

In accordance with what has been previously stated, a more strict 
consideration of the problem, appears desirable, that is, an explanation 


*"K Teorii Rasprostranenia Plameni." Zhurnal Fizicheskoi Khimii 
(USSR), T. 22, 1948, pp. 27-49. 
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of those general conditions that must be satisfied by an explosive 
mixture for the possible existence of the phenomenon of flame propa- 
gationj determination of a general method of computing the flame 
velocity for an arbitrary relation among rate of reaction, tempera- 
tiire, and concentration; and finally, determination and evalua- 
tion of the errors of the aforementioned approximate method (refer- 
ence 3 0). 

The present paper bears greater resemblance in its form to a 
mathematical than to a physicochemical investigation. The subject, 
however, is actually the physicochemical problem of the conditions 
under which the chemical reaction proceeds dtiring combustion. The 
theory developed permits studying the rate of the reactions occurring 
at temperatures of 1500 to 3000°. C for a time of the order of a 
millisecond to enable the measirrement of the flame velocity. The 
mathematical apparatus provided is not complicated; the computation 
work is simularly not extensive particularly if it is compared with 
the difficulty of the direct experimental realization of the indepen- 
dent combustion of the gas and the direct measurement of the rapid 
homogeneous reaction under conditions approximating the flame 
conditions. 

At the present time, investigations are being conducted and 
prepared for publication that concern concrete chemical systems that 
will extend our knowledge of the kinetics of the most important 
reactions to the region of high temperatures; and will provide a 
clear picture of the process of combustion of such propagated explo- 
sive mixtures as carbon monoxide - oxygen - nitrogen. Accordingly, 
a review of the present state of the theory that presents not only 
the final results but the entire basis of the theory, emphasizing 
the fundamental physical problems of the existence and properties 
of the solution in the general case, would seem desirable. 


1. Equations of Heat Conductivity and Diffusion in a Flame 

A two-dimensional flame is considered and the eqijation of heat 
conductivity giving the distribution of the temperatures and the 
equation of diffusion giving the distribution of the concentrations 
are set up; the total number of equations of diffusion is equal to 
the number of reacting substances and the products of the reaction. 
Only one equation of diffusion is ivritten for the reacting substance 
of which there is a deficiency in the explosive mixture. In con- 
trast to the classical problems of heat conductivity and diffusion, 
the liberation of heat and the consumption of the reacting substances 
that depends on the process of the chemical reaction are considered 
in the equations . 


NACA TM 1282 


3 


The following symbols are used; 
t time, (sec) 

X coordinate, x-axis perpendicular to plane of flame, (cm) 

p density of mixture, (g/cm^) 

u velocity of motion along x-axis, (cra/sec) 

T temperature, {deg) 

c specific heat at constant pressvire, (cal/(deg) (g) ) 

k heat conductivity, ( cal/ (cm) (sec) (deg) ) 

a relative concentration of reacting substance in grams of 
substance per one gram of mixture (nondimensional) 

2 

D coefficient of diffusion, (cm^/sec) 

velocity of chemical reaction, (g/(sec) (cm^) ) 
h thermal effect of reaction, (cal/g) 

X constant of thermal diffusion (nondimensional) 

The equation of heat conductivity has the usual form 


5t S 


3t St , ^ 

- wc + hi 


( 1 ) 


The equation of diffusion may be written 


p ^ = A. 

St Sx 


D 


Dp ^ ^ Dpa . ^ 

Sx T Sx 


] - 


up 


Sa 

Sx 


- <1» 


( 2 ) 


With change in temperature, the density of the gas is found to 
change greatly. The composition of the mixture is characterized 
by the magnitude a. The absolute value of the concentration, equal 
to ap (g/cm^), changes because of gas expansion, diffusion, and 
reaction. 


In equation (2) only processes that change the composition 
are considered and the expansion is excluded. The magnitude Dp (Sa/Sx) 
represents the mass flow of the substance carried along by the dif- 
fusion; in an Isothermal diffusion, the flow is generally confeidered 
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to be proportional to the gradient of the absolute concentration 
pa, that is, in the case D [§(pa)/d^ . In this case, however, 

the dei^sity p is practiQally constant and whether the mass flow 
is written as D [S(pa)/S^ or Dp(Sa/Sx) is immaterial. In 
the nonisothermal case, the second form must be definitely chosen 
because, in the absence of thermodiffusion (separately considered 
by the second term •)r(Dpa/T) (5 t/Sx) ) in a mixture of constant 
composition, there is no diffusive flowj that is, the flow is equal 
to 0 for a = constant, even if pa ^ constant in this case. 

Equation (2) may be obtained by a transformation of the equation 
for the absolute concentration 

A (pa) = A ('dp^+T^-^VuA (pa) - pa A - $ (3) 

Sx 5x ^ Sx T Sx / Sx hx 

by considering the equation of continuity 

5p ^ 5p p — 

5t Bx Sx 

The magnitude 4> represents the rate of consumption of the 
substance in the process of the chemical reaction; the rate of 
heat liberation in unit volume is equal to h<J). $ is a f\inction 

$(a,T) where = 0 for a = 0. 

By considering the propagation of the flame, the velocity 
of motion of cold gas Uq is sought for which the flame front is 
stationary and the entire distribution of T and a is independent 
of time. It is evident that Uq is equal to the velocity with 
which the flame is propagated over the stationary cold gas but it 
is more convenient to consider the flow and the distributions of 
T and a as being independent of time. As the temperature rises 
in a gas flow, p drops and u increases in such a manner that 
the product pu is constant at any point and equal to Pq^O' 

In equations (l) and (2), the left side (S/St) is equal to 0; 
the symbol S/Sx is replaced by d/dx. The following boundary 
conditions are imposed on T(x) and a(x) : 


for x=-® T=Tq a= 
for X = + 00 a = 0 



(4) 
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(because for a 0 and ^ ^ 0 , the condition a = constant 0 does 
not satisfy the equation) . 

It is likewise evident that for x = ±oo 


dT da 
dx dx 


( 5 ') 


2. General Properties of Temperature Distribution 


The general properties 
respect to x from, -oo to +c» 


of T(x) are found by integrating with 
in both equations. 


— k — 
dx dx 


upc 


dx 


+ 


h4> 


0 


( 6 ) 


By setting 


_i -up^-*=o 

dx \ dx T dx / dx 

T(+t») = Tg; 


upc (Tg - Tq) 


-^±(X 

h $dx 

\»/ — OO 


( 7 ) 


( 8 ) 


upa^ = 


-f“OD 


$dx 


har 

^2 = To 


( 9 ) 

^10) 


Equation (lO) agrees 'vri.th the result of the elementary computation 
of the combustion temperature by the energy balance. By considering 
equation (6), it is show that T cannot pass through a maximum 
higher than Tg because the curve passing through the 

maximum, however, dT/dx < 0 and (d/dx) (kdT/dx) >0, which according 
to equation (6) requires <I> < 0. 


In reference 1, incorrect assumptions are encountered, that is, 
for a rapid reaction of the gas heated to a certain temperature 
T-g higher than Tq, a sudden increase in temperature occurs and 


T 


In the case of variable 


c, set 


H = 



cdT. 
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the following temperature is attained; 




+ 


aoh 


c 


3 » 111 

higher than T2 = Tq Q « Actually, if $ is very large in a 

very narrow layer so that ^ is finite, then d^T/dx^ in this 

layer is very large, that is, in the limit as $-»oo, the magnitude 
dT/dx undergoes a discontinuity hut T is continuous; an angle appears 
on the temperature -distribution curve (fig.l) . It is also evident 
that the converse is really true; namely, that at all points at 
which $ is finite or equal to 0, not only T but also dT/dx 
cannot undergo a discontinuity. 


3. Relation Between Concentration and Temperature in Mixtiore of 
Gases of Approximately Same Molecular Weights 

A mixture of gases of nearly the same molecular weights is next 
considered. In this case, it follows from the kinetic theory that 
practically no thermodiffusion occrrrs (y = O) and that the dif- 
fusion coefficient is related to heat conductivity by the expression^ 


k = cpD (11) 

The variable H is introduced (heat content or enthalpy) , uniquely 
related to T as follows ; 



dH = cdT 



( 12 ) 


^At room temperature, the heat conductivity of oxygen is 
5.9 X 10"^, nitrogen is 5.2 X 10"^, carbon monoxide is 5.1 X 10“^, 
and hydrogen is 40.7 X 10"^. The product cpD for diffusion of 
oxygen - nitrogen is 5.8 x 10"^, hydrogen - oxygen is 23 x IQ-^, 
hydrogen - nitrogen is 21 x 10“5, and hydrogen - hydrogen is 
38 X 10“^. 
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The equations of heat conductivity and diffusion assume the similar 
forms 

_d_lsJS-up^ + h$=0 (13) 

dx c dx dx 


_iiS^_up^-$ = 0 (14) 

dx c dx dx 

from which 

— — — (H + ah) - up — (H + ah ) = 0 (l5) 

dx c dx dx 

thus, by considering the boimdary conditions (S') that will give 


_d_ 

dx 


(H + ah) = 0 for x = + °° 


(15’) 


the unique solution 


H + ah = constant = Hq + aQh 


(16) 


X = - OB 

H = Hq 
a = Uq 

X = + Oc 

a =t 0 




( 17 ) 


H = H2 = Hq + agh J 
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The concentration is found to be uniquely related to the temperature 
(in the case of constant specific heat linearly: cT + ah = constant 
= cTq + aQh) in such a manner that the total energy^ of a unit mass; 

the sum of the thermal and chemical energies is constant. 

Levis and von Elbe (reference 9) were the first to postulate 
the constancy of the total energy in all the intermediate layers in 
the flame by the assumption of a chain reaction that is brought 
about by the diffusion of active centers. From the previous dis- 
cussion^ it is evident that the constancy of the energy depends 
not on the mechanism of the reaction but on the ratio between the 
diffusion of the fuel and the heat conductivity of the mixture. 

As will be shown, the total energy is not constant for Dcp ^ k 
in the region of combustion, but has a maximum or minimimi value. 

In considering the case Dcp = k, the velocity of the reaction 
4> (a,T) may be expressed for known Uq, Tq as a function of one 
variable magnitude a or T, or H because of the relation among 
the composition of the mixture, the temperature, and H. 

The problem thus reduces to the consideration of the one equation 


_d_kjiH.up.!iH + h$(H) = 0 (18) 

dx c dx dx 

in which *(h) is obtained from 'l>(a,T) by the substitution of 
T and a from equations (12) and (16) so that for H = H 2 ^ $ = 0 

because for H = a = 0. Equation (18) may be solved imder 
boundary conditions (l7) and 


X = - t» 

H= Ho 

X = + Oo 


H = E 






2 J 


(18a) 


%ore accurately, the enthalpy (heat content) because the com- 
bustion occurs at constant pressure . See definition of C. 
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The product up is a constant magnitude that cannot, however, be 
considered as knoim; it is necessary to determine that value of up 
for which the equation has a solution satisfying the boundary 
conditions . 

4. Condition of Existence and Uniqueness of Solution 


The following nondlmensional variables and parameters are 
introduced^ : 


^2 - H 

H2 - Ho 



where K^, are constants independent of a, H, and other 

magnitudes varying in the flame front. All finite results of the 
computations, particularly the magnitude of the velocity u, do 
not actually depend on the choice of K^, and Cj^. Because 

1 _1 
^~^m ^ found farom equation (20) that m~^*jjj ^ and that 

1 

the velocity u~ does not depend on ^ hence, the choice 


^This convenient system was proposed by S.M. Painberg in con- 
nection with work on the combustion of ozone. 
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of and is arbitrary. For example^ set 

Kj^ = k (a = 0, T = Tg) 

Cj^ = C (a = 0, T = Tg) 

^m= 



(19a) 


It is noted that for Dcp = k. 


z 


a 



By substituting and noting that 


A = ^ A 

dx dx dH 


the following equation is obtained in place of equation (l8): 


y g + my _ cp(z) = 0 


( 20 ) 


with the boundary conditions 


z = 0 
y = 0 
z = 1 
y = 0 




> 


J 


( 21 ) 


The conditions of the existence of a solution are considered. In 
the plane z^y the line y 2 ^ = ip(z)/m is drawn on which dy/dz = 0 

(fig. 2)^. For y< y^^ dy/dz > 0. From this condition it is evident 


^y^^ is not the integral ciirve of equation (20) but only one of 
the isoclinal equations. 


t\ f' 
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that if cp(z) >0 in the entire interval 1, the curve 

satisfying the condition y « 0, z 0 lies everywhere above the 
axis of abscissas and never (for any m) falls on the point 
z = 1^ y - 0. This condition exists because near the point 
z=*l, y« 0, dy/dz > 0 and there are no curves from the region 
z < 1, y > 0 that pass throiagh this point. Hence, in order that 
a solution exist, it is necessary that Cp becomes 0 for z$ 1. 

The case considered first is where an even more strict condi- 
tion is satisfied. 


cp(z) a 0 for z > p, where 8<1 (22) 

(Note fig. 3.) Also, let 9(z) > 0 for z < p. The existence 
and uniqueness of the solution will be proven. The general theorems 
on the tiniqueness of solution of a given differential equation are 
insufficient because the value of the coefficient m must be found. 

Consider y (z,m), a function of two variables, as satisfying 
constant ra of equation (20) and the boimdary condition y (0,m) = 0. 
The following equation is set up; 


Sy 

By dividing equation (20) by y and differentiating with respect 
to m 


(23) 




m 


cp (z) 




dz [y (z, m^ 
z = 0 


2 m 


y^ = - 1 


(24) 


y = 0 

In the region 0$.z$p, cp>0 and y> 0; it follows that y^^^ < 0. 
Set 


-ym (3^m) a e 

In the region P ^ z < 1, dy^/dz » - 1 so that 


(25) 


- yja (I, m) = € + (1 -P ) > 1 - P 


(26) 
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Hence, dy/dm for z = 1 is alvays negative and an absolute value 
greater than a certain finite magnitude. 

For m = 0, the initial equation can readily be solved as follows: 



As has been shown, as m increases, y (z = l,m) decreases mono- 
tonically and its derivative with respect to m is always finite. 

Thus, there is always one and only one itiq for which 

y (z = 1, ruo) = 0 (29) 

This value for which both conditions (y = 0, z = 0 and 

y = 0, z = l) are satisfied will give, by formulas (l9), the required 
flame velocity. 

The practical method of finding m for a given <?(z) and a 
given arbitrary value consists of finding the corresponding 

y (z = 1, m-^^) by numerical integration of equation (20) from the 
point z=0, y=0. If y(z = l, m^) ^ 0, another mg, mg > m^^, 
is chosen and the nimierical integration is again carried out^ if 
y (z = 1, m^) <0, m^^ is chosen. Several values of y (z = l) 

are found for various m = m^, mg, m3 and are chosen so that some of 


6if the other sign of the root is taken, the sign of m changes, 
that is, the propagation is considered in the reverse direction but 
all the physical conclusions remain valid. 
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the values found for y (z *= l) are positive and some negative. By 
laying off on a graph y (z = l#ni) as a function of m, the value 
of niQ for which y (z = 1, ioq) = 0 is found by interpolation. 

As a result, the velocity of coinbustion is obtained by the formula 



where represents a nondimensional number determined by the 

fimction cp , that is, by the cirrve of $k/c as a function of H 
(or of the temperature) . The dependence of the flame velocity on the 
absolute values of the reaction velocity and the heat conductivity 
is given by the root of formula (30) . In considering the dependence 
on the composition of the mixture, it is necessary to take into 
accoiint that aQ is an important factor in the expressions for 
$ and because in passing from 0(a,T) to ^(H), Sq entered 

formula (16) relating a and H. 

By knowing hiq and the corresponding integral curve y (z,mQ), 
the law of distribution of H, T, and a in space as a. function of 
X can also be found. For this determination, a single quadrature 
is sufficient that can always be carried out numerically 



where the subscript m denotes all the scalar magnitudes (see equa- 
tions (l9a)). The expression \mder the radical represents a constant 
magnitude having the dimension of a length and the expression in 
parentheses is a nondimensional ratio depending on z. 



X 


(32) 
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By the given condition, x = - <*> for z ** 1 and for x = +» , 
z = 0. Near z = 1, y (z) = m (l - z) that gives, after substi- 
tuting the value m, the known Michelson solution 


X = — In (T - T„) + constant 


/V \ 

l./^ T UC1.X1 I ' I 

° \<=0<=oV 


T = + constant ’ exp 


(33) 


Thus T Tq only asymptotically as x - oo. The behavior of the 
solution near z = 0, that is, where the reaction ends, depends on 
the form of y (z) near z = 0, which in turn is determined by the 
expression of Cp (z) for small z. 

If <p(z) increases near z = 0, according to the law 


cp(z) = constant z’^ 
n >1 


it is easy to show that 


(34) 


y (z) 


m 


so that if the Integral cp(z) diverges, the integral 

/ dz/y(z) likewise diverges. In the absence of diffusion and 

w 0 

heat conduction for an adiabatic reaction in a closed vessel, 
the magnitude z varies according to the equation 


(35) 


H = - A q)(z) (36) 

at 

where A = (l/a^) ($k/c)jj^ (c/kp) in the limit for small z approaches 
a constant value. Thus, if the Integral / dz/cp(z) diverges, 

then, under adiabatic conditions, the chemical reaction comes to 
an end only asymptotically at t -► oo. Correspondingly, in consi- 
dering the problem of flame propagation, the end of the reaction 
a = 0, T = Tg, and H = Hg is reached only in the limit at an 
infinite distance from the flame front. 
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The reverse case for small z is considered. 

n - 1 

y (z) » constant • z ^ 


n + 1 


< 1 


(37) 


In this case. 


cp(z) 


constant 


n < 1 


n 


Thus, in the case where the integral 


(38) 


integral 


dz/y(z) 


dz/ cp(z) 


is finite, the 


is likewise finite and converges nesu* z = 0. 


Physically, the chemical reaction concerned goes tg completion 
after a finite time under adiabatic conditions (in a closed vessel). 
In the presence of diffusion and heat conduction, the chemical 
reaction in the flame is likewise found to go to completion at a 
finite distance from the flame front (more specifically, at a finite 
(^stance from the place of most intense reaction) . By computing 

for that value 


dz/y(z) (for example, arbitrarily giving 
^ for which <f>(z) = maximum), Xj^ 


X = 0 
is found for which 


X = Xj, 

a = 0 

H = H. 
c 

T = T, 


dx 




( 39 ) 


J 


The distribution of the concentrations and temperatures for x < Xj^ 
is given by the integral (32) . For any 
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X ^ 

a = 0 



> 


(40) 


Thus, in the case of finite 



J 

a singular point at 


X 




appears in the distribution of a(x), T(x), and H(x). At this singular 
point, two different solutions, relations (32) and (40), meet and 
the values of the magnitudes T, a, and H and their first deriv- 
atives with respect to x are the same on both sides'^. The second 
or even higher derivative, however, may suffer a discontinuity, 
particularly if 


qp(z) = constant 
- 1 < n •< 1 


near z = 0, 

2 

Tg - T = constant (Hg - H) = constant • a = constant (x^ - x)^ “ 

(41) 

where from the conditions® imposed on the exponent, it follows that 
the degree of (xj^ - x) is positive and greater than 1 so that 
Tg - T = 0 and dT/dx = 0 for x -*• Xj^. 


"^The first derivatives can have a discontinuity only for 

®The condition n < 1 follows from the requirement of the 
finiteness of P dz/cpj otherwise, as has been shown, there is no 

Jo 

singular point on the curve. The condition 1< n follows from 
the requirement of the finiteness of J cpdz. As will be shown, 

Jo 

if this integral diverges, the magnitude m and the flame velocity 
are infinite; however, n< 0 already leads to infinite velocity 
of reaction for a-» 0, which is physically unreal. 
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The relations (32) and (40) thus give a complete solution sat- 
isfying equation (l8) and ho-undary conditions (l8a) . The fact that 
the "boundary condition H - Hg is satisfied not only for x = » 
hut also for all x > xj^ does not contradict or prevent the solu- 
tion from being considered correct. The part played hy J^dz/ q>(z) 

was pointed out to the author hy "7. I. Skohelhin hut the author cannot 
agree with the opinion of Skohelkin that for the existence of a 
solution it is necessary that this integral diverge. 

A solution has heen shovm to exist only in the case where for 
z = 1, cp = 0, that is, in the case where the reaction velocity is 0 
in the initial mixture; for T = Tq, a = a^. The solution has heen 
found mder the still more rigorous condition ^22). This result 
is in no way an indication of the insufficient generality of the 
theory and the inapplicahility of the theory to any case of flame 
propagation. On the contrary, this result entirely corresponds to 
those physical conditions that must exist in order that a strictly 
stationary process of flame propagation can occur. Not in every 
system in which the process of a chemical reaction is thermodynamically 
possible will flame propagation he observed, that is, the propaga- 
tion in the space of the region of chemical reaction. For this flame 
propagation, it is necessary first of all that the heat given out 
during the reaction accelerate the process of the reaction in the 
neighboring layer. This condition, however, is still not suffi- 
cient; it is also necessary that in the Initial state, that is, in 
the state in which the mixture is found at x = - at an infinite 
distance from the flame front at the Initial temperature Tq, the 
reaction velocity is strictly equal to 0. If this condition is 
not satisfied, a certain finite time may he found during which the 
mixture, which at the initial instant is in the initial state, reacts 
in the absence of external actions. This phenomenon can he illus- 
trated hy Imagining an explosive mixtirre in an infinitely long 
pipe. If at the Initial instant the mixture throughout is in the 
initial state except at the origin of coordinates where the flame 
is situated, the time required for the flame to reach a definite 
point in the pipe is equal to the distance of this point from the 
origin of coordinates divided hy the flame velocity. An element of 
volume of the mixture with a finite reaction velocity at Tq and 
at a sufficiently long distance from the origin reacts before the 
flame reaches it. As a result of the reaction as the flame is 
propagated, the properties of the mixture situated ahead of the 
flame front at Infinity vary. Hence, for $(TQ,aQ) ^ 0, the 
required regime of steady flame propagation is actually physically 
impossible and the absence of a solution of the equation in this 
case entirely agrees with the physical picture. 
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The peculiar mathematical difficulties aj-ising in the consi- 
deration of reactions for which Cp(z = l) becomes 0 , for example, 
pixsportional to (l - z) , is considered. In the case of heat propa- 
gation, there is no significance in investigating the flame velocity- 
on the assumption that for T near Tq, the reaction velocity is 
proportional to (T - Tq) because it is experimentally possible to 
vary Tq so that the flame propagation does not change. If for 
some Tq = IqI' ^ egual to constant (T - then by this 

assumption, the case previously considered, Cp(z> g) = 0 where 
g < 1, is obtained for Tq 2 < Tq^> for Tq 3 > Tq^, Cp(l) / 0 is 
obtained, that is, the impossibility of a steady propagation. Tlius, 
in the heat propagation the assumption $ ~ (T - Tq), $ ~ (l - z) 
is a random one and unsuitable. 

From general molecular-kinetic considerations it follows that 
at a low temperature near Tq the velocity of the chemical reaction 
is finite and different from 0 but very small. In this sense, 
the existence of a steady flame propagation is an approximation 
because $ (uq, Tq) / 0, but a very good approximation because 
^)(aQ, Tq) is very small. By the law of Arrhenius, the velocity 
of the reaction at room temperature is of the order of e"^® or e"^*^ 
(for -unity, the velocity at 3000°C was taken) . 

The curves of (p/m and y as functions of z are shown 
schematically in figure 4. For z = 1, cp / 0 but for z > 

(p/m4C y. From a consideration of the equation, it follows that to 
obtain a solution it is necessary to set 'p(z) = 0 for z >g; 
it is also easy to show that all numerical results (ra, u, and the 
curve T (x)) depend very little on the choice of g (fig. 4^ 
or the larger P 2 ) ‘PCz)'^ 1 for z> 

Thus, although the operation itself of cutting off 9(z) at 
z > g contains an arbitrary element in the choice of 3, the cut-off 
is required and arbitrariness in the choice does not lead to arbi- 
trariness in the results. In this sense, the theory may be considered 
satisfactory. 


5. Isothermal Chain of Flame Propagation 

The cut-off by means of multiplication by (l - z) in the heat 
propagation of the flame is, as has been shown, a random one. In 
considering chain and autocatalytic reactions, however, the assumption 
that the velocity of the chemical reaction is proportional to the 
concentration of the final product formed in the reaction is of 
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interest. For example, if the transformation of A and B (their 
concentrations are denoted by a and b) is assimed to proceed with 
the collisions of the molecules of A and B according to the 
equation 

A + B = 2B 

then, by assimming a constant temperature, the possibility of flame 
propagation is obtained if the substance B transformed during the 
reaction diffuses in the neighboring layers where the pure substance 
A is contained and results in a chemical reaction of the two sub- 
stances, in which the velocity constant is denoted by C. An equa- 
tion identical to equation (20) is obtained with 


(t(a, b) = C • a • b ^ 


*m = Cao2 

^(z) = z (l - z) J 

In connection with the problem of the propagation of a chemical 
reaction, this problem was considered in the unpublished works of 
B. N. Skalov and 0. M. Todes at the Institute of Chemical Physics 
during 1937 and 1938. The results of this work are given in refer- 
ence 12. Attention has also been called to the very interesting 
work of A. N. Kolmogorov, I. G. Petrovskii, and N. S. Piskunov 
(reference 13), which was published in 1937. In this work concern- 
ing the equivalent biological problem, all the results of this 
problem obtained later at the Institute of Chemical Physics are 
presented in an accurate and clear form. 

It has been shown that equation (20) for cp(z) given in 
equation (42) has a solution satisfying the boundary conditions 
for all 

®( 43 ) 


^Solutions are also possible -for m < 2 but these solutions 
require the existence of the region z > 1, a > a^, and b < 0 in 
which $ < 0, thus indicating their complete physical unreality. 



> 


(42) 
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The physical meaning of these solutions has been clarified. For 
cp(z), given by formula (42), even in the absence of diffusion) 
the initial distribution a and b in space may be given satis- 
fying the condition a = a^ and b = 0 for x-» which gives 

the propagation of the reaction in space. As quoted from reference 13; 
"The apparent displacement of the substance from the left toward the 
right will actually be brought about here by the Increase of its 
density at each point, entirely independent of the course of the 
process at other points." The time of the reaction according to the 
equation 

^ = Cab = a (ao - a) (42a) 

depends logarithmically on the initial value s.Q-a. that approaches 
infinity for s,q~ 8- -> 0 and a-> a^. Given the distribution 


a = Uq (l - e’^) (for x < O) 


(44) 


that satisfies the condition a -* Sq, x a linear dependence 

of the reaction time on the coordinate is obtained, which leads to 
an apparent displacement with constant velocity. The smaller n 
is the greater is this velocity. By considering the nonstationary 
problem that for the particular case (42) has the form 


da _ 

dt 


D ^ 5 a (a^, 

dx2 P ° 


a) 


(45) 


Kolmogorov, Petrovskii, and Piskunov showed that as the time t->®, 
the portion of the curve a(x, t) on which the initial part of the 
change of a from aQ to 0 occurs is displaced in space with a 
velocity approaching the minimum value that satisfies the steady- 
state equation 


u = 



(45a) 


and the distribution a(x, t) approaches the solution of the 
steady-state equation. The physical result is confirmed that in 
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any case of local inflammation the reaction will be displaced 
with a velocity according to equation (45a) and that larger velo- 
cities are physically unreal. ^ 

According to Skalov, the limiting maximum value of the velo- 
city in the nonsteady problem may be found by elementary methods by 
considering the transformed equation 


Sa _ -n 5^a 


St 


Sx^ 


- C a (a _ a) 


(46) 




that agrees with equation (45) for 
sidered throughout for the initial conditions 
for a small region near x = 0, where a aQ. 

C 

p ^0^ 

a = Uq - qe 


Equation (46) is con- 


t = 0, a = 
By setting 


except 


(47) 


the following equation of diffusion without sources is obtained for q: 


St " Sx2 


a = A • 




> 


j 


(48) 


a 


a 


0 


^ _1_ eP - 4Dt 


(49) 


The relation between x and t is found at the point at which 
a = a' < aQ, that is, the law of displacement in space for the given 
concentration a' 


C 

P 


a^t 


4Dt 


In (bq - a ' ) 


In A 


+ i Int 


(50) 


^^The unreality of a greater u may also be established by the 
method of a cut by considering the following relations: $ = 0, 

aQ-a <e, $ = Ca (aQ-a), and aQ-a> g and by finding in this case 
the only value of m and the corresponding value of u. For €-*0, 
u 2 y^(c/p) Da^. 
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In the limit as t 


x2 = 


X = t 


4DCaQt2 

P 



(51) 


which was required to be shown. The substitution of (c/p) a(aQ- a) 
by (c/p) aQ(aQ- a) is an overevaluation and can give only an 
increased but not a decreased velocity. 


6. Approximate Fomiilas for Reactions Strongly 
Dependent on Temperature 

Considering the theory of the thermal propagation of a flame, 
approximate expressions of the flame velocity are considered for 
the case where the velocity of the reaction increases rapidly with 
increasing temperature so that the function $ of equation (18) 
(in which a and T are already expressed in terms of H) has a 
sharp maximum for H near Hg. 


In equation (20), ^(z) is shown to have a sharp maximum for 

z <s 1. This case is of most interest for typical chemical reactions 
of combustion. 


According to assvmiption (22) that 9(z) is different from 0 
only for z < p, a strict evaluation can be given of the upper and 
lower limits of the magnitude m in such form that for 6“^0, they 
approach the same value that is the limiting expression for m for 
small P . 


The entire range of integration of z from 0 to 
divided into two parts: from 0 to p and from p to 
second part 


1 is 

1 . In the 


<;p = 0 


dz 


- m 


y = y(P) - ra(z - 3) 




y(l) = 0 = y(p) - m(l - 3) J 


( 52 ) 
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hence, it follows that 


m 


y(P) 

1 - P 


(53) 


so that finding m reduces to determining y(P), where an exag- 
gerated value of y(8) will give an exaggerated value of m. In 
the first range, the exaggerated value y ^ is obtained, rejecting 

IU3«X 

the negative term my: 



my -c y 


max 


^ymax 

dz 


Cp 


(54) 



where u is an auxiliary variable of integration and I denotes 


the integral 



The limits of integration may also be extended 


to luiity without markedly changing the integral because the case of 
the function cp is considered as having a sharp maximum at small 
z < p . 


Thus, by taking 



(56) 


the upper limit of m is obtained 


%ax = 


-> (1 + P + . . . ) for PCI 


(57) 
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In order to find the lower limit y . of y, the equation 

min 

is Witten in the form solved for dy/dz and the dominating value 
^max from equation (55) is substituted on the right side for y. 


tiy m 

^ - m > 

y ymax 


dy . 

CD '^min 

-t — - m = 


dz 


y >y. 


• =r 

min J 


q?(u) du 


0 


'max 



- mz = . / 2 


J^Cpdcp - 


mz 


(58) 


'min 



Thus; for small p, the upper and lower limits coincide and 
give the limiting expression 



(59) 


Chemical kinetics generally lead to the expression 


$(a, T) = constant • a^e"-^/^'^ 


(60) 


2040 
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which after transformation gives 


(p(z) = Cz^e ^ 




"^0 

x = 1 <1 


J 


( 61 ) 


This expression may be approximately represented by making use 
of the assiunption of Frank-Kamenetskii (reference 17) in the form 


cp ( z) = Ce 


n 

z 


e-xXz 


(62) 


which in agreement with equation (59) (considering the fact that 
e”^^<Cl), gives 



(63) 


Expressions (59) and (63) represent those approximate solutions 
that are given in references 10 and 15 v/here the physicochemical 
consequences of these results are given in detail. The less accurate 
but clearer derivation of these formulas directly from equations (6) 
or (18) in the coordinates x,T or x,H will not be repeated; as 
has been shown, a strict mathematical analysis confirms the correct- 
ness of this result. Expressions (60) and (62) do not become 0 
for T = Tq and z = 1. The nimerical computation shows, however, 

that they are then very small (as exponents of large negative 
numbers). The principal difficulties connected with this and the 
method of overcoming them were previously explained. Evidently in 
cutting off 9 (z) in the region where cp(z)4C 1, the part played by 
the magnitude 6 ( 9=0 for z > p) will be played by the magnitude 
where |J. is of the order of several units. 

In this manner, expression (63) and corresponding formulas of 
the previous work are correct in the limit when x\»i (which 
corresponds to the case p <3^ l) , that is, for 

A (Tg - Tq) 


»1 


(64) 
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as has heen shown in references 10 and 15. 

By the method of successive approximations, a correction of 
the order l/^X., l/^ and their powers can also he obtained in 
expression (63) (see reference 16 for an example by P. G. Smirnov). 
In view of our present knowledge of the velocity of a chemical 
reaction and a number of idealized assumptions (for example, the 
absence of heat losses) , however, such exactness is of no practical 
interest. 


7. Case of Mixture of Gases of Different Molecular Weights 

Those cases where the identical relation (16) between the con- 
centrations and the tfemperature does not hold will be briefly 
considered. 

Such a relation does not hold if: 

a) Dcp ^ k and y ^ 0, which holds true in a mixtvire of gases 
differing considerably in their molecular weights. 

b) Several (parallel or successive) chemical reactions take 
place in the flame and the relations between their equations and 
the velocity are such that the concentrations of the different 
chemical substances are not connected by algebraic relations. 

In the first case, equations (6) and (7) must be considered, 

A nondiraensional transformation analogous to equation (19) is 
again introduced. For exactness, '^ 2 ) taken for the 

unit of reaction velocity; K = k(0,T„) and C = c(0,T„)the values 

m u in o 

of the combustion products, are taken for the \mits of heat conduc- 
tivity and heat capacity. 

The nondiraensional relative concentration a and the other 
nondiraensional magnitudes already employed are also introduced. 
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a = 




z = 


Hg - H 






(p(a,z) = 


"2 “0 
k (aT) $ (a,T) C 


m 


C (a,T) K $ 

' ' ' mm 


> 


^ “ k dx 


m = pu 


c dz /C^ (Hg _ Hq) 


K i 

m m 


(Hg - Hq) 




(65) 


H 


2 


H 


ha^ 

0 


Expressed in these variables, the heat conductivity and- diffusion 
are transformed into 


y 


dz 


+ ray - 


cp (a,z) 


0 


( 66 ) 


where 



da 

dz 



day + 



cp (a, z) = 0 




^ (a,z) 


Dcp 

k 


9=9 (a,z) 


rDp (Hg - Hq) 
kT 


(67) 


The boundary conditions are z = 0, y = 0, a = 0, z = 1, y = 0, and 
a = 1. The nondimensional magnitudes and 9 characterize the 
ratio of the diffusion to the temperature conductivity and the 
thermodiffusion, respectively. In general, both magnitudes depend on 
the composition and the temperature, that is, they depend on the 
nondimensional variables a and z but they have neither singu- 
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laritles nor zeros in the entire range of variation of a and z. 

At all times , d>0) iJ>l for a mixtxire containing a small amount 
of a light gas (for example, Hg + Og or Eg + Erg in an insufficient 

quantity of hydrogen) . As a rule, in mixtures in which ^ > 1, 
y < 0, and 0 < 0 (in a thermodiffusion equilibrium, a light gas 
generally concentrates in the hot part of the vessel) . 


For ^ = 1 and 0=0, the system has the obvious solution 
a = z, which \ras previously used. In the general case of arbitrary 
0, and a method may be proposed analogous to that previously 

described that consists of the numerical integration of the equations 
for different m at the initial stream conditions for z = 0 and 
the choice of interpolation of such a value of m for which the 
solution satisfies the conditions for z = 1. 


The difficulty is that according to the general theory for the 
integration of an equation of the second order for a, it is neces- 
sary for z = 0 to give, in addition to a = 0, the derivative 
da/dz. An explanation is also necessary as to whether for one m 
both conditions a = 1 and y = 0 will similLtaneously be satisfied. 


By dividing equations (66) and (67) by y and integrating, the 
following equation is obtained for the solutions satisfying the con- 
dition a = 0, y = 0 for z = 0: 


i 

r* n icL da,| \ 

y(l) + m = r - dz = ^ - iSy - 0a(l) y(l) + ma(l) 

0 


( 68 ) 


Hence, that value of m for which y(l) * 0 also sat^|fies 
a(l) = 1 if y da/dz = 0 for z = 0 and for z = 1. 


^^On cp are Imposed the conditions <p = 0 for a=0, i?)=0 in 
the neighborhood of a = 1, z = 1, that is, $ = 0 for T, a near 
Tq, Uq in the initial mixture as is required for the existence of the 

propagation regime (see preceding discussion). 

1 ? 

The satisfying of y = 0 for z = 0 and z = 1 is not suf- 
ficient for this because -da/dz may approach infinity. 


Ill 
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Near z = 1, a = 1, and y = 0 
the equations reduce to the system 

dz 


y = m (l - z) 




Thus, hy whatever law the curve a(z) approaches the point z = 
for that value of m at which point y = 0, 


d X da 

dz ~ dz " dz 


^ 9a (1 - z) f 


By retaining the principal terms (for small 1 - 

1 


d 

dz 


I 


_ ^ 

z dz 


e 


1 - z 
1 


da 

dz 


e 


Y + constant (l - z) 




- 1 


7 


in the region in which q> = 


d.a \ da _ „ 

y _ = ^(1 . = 0 

for z = 1 also in the case (^ < l) where da/dz -> <» for z = 1 
For obtaining m, it is therefore sufficient to give 


z = 0 




y = 0 
a = 0 



y 


da 

dz 



^^The values of ^ and 0 in equation (70) are taken for 
a = 1. 


0 , 


(69) 


(70) 


1. 


(71) 


z ** 1, 
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and find the value of m for which the solution satisfies y =? 0 
for z .= 1; the conditions a » 1 and z = 1 will he automatically 
satisfied. 

By assuming qp for small a in the form cp = z) where 

\(r(0) ^ 0, the following relations are obtained for small z« 1: 


a = z for n > 1 


a = £ for n < 1 


(72) 


-1 


a = z 




m2 

~2Ho) 


form: 


for n = 1 

As a rule, in the combustion of gases, the function has the 


cp(a, z) = const = const a’^e 


-A\z 


(73) 


where A * E/RTp > 1. 


The coefficient A usually lies between 7 and 20 and 
\ t= (Tg - Tq)/T 2 , 1 - \ = ^q/'^ 2 geiierally lies between l/4 and l/lO. 

For z = 0 and qj = 0, the following relation is obtained for very 
small z : 

z < -^ « 1 

AX 

Cp increases with increasing z because of increasing a but with 
further increase in z,<P rapidly (exponentially) drops. The region 
of maximum Cp near z = n/AX<C 1 is of principal effect for the 

integral J^qpdz. 

From estimates (57) and (58) previously made, it follows that 
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i. 


i 


I i 


In the principal region near maximum 9 and maximum y 

my3f2 y ^ ( 74 ) 

whence, in the region of chemical reaction, the principal terms of 
equation (67) are 


qp(a. 


z) 



da 

dz 


(the terms rejected are less than the remaining terms in the ratio 
n/A\) therefore, it follows that in this region 


d o da dx 

y — Y — 

dz dz dz 


dz " 


> 


(75) 


IX 'St — I 

Independent of the value of m determining the ratio of a and z 
in the limit for very small z. 


The result given, valid in the limit for very large heat of 
activation (greater than A), was pointed out by L. L. Landay. 


For this limiting case, the following simple formula is obtained 
for the flame velocity: 


where 


m = 


_L 

r 


qpdz 


(76) 


2 

cp = 9 (a, z) = 9(0- = z) 


HIM ■mil II 
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Thus, in the case of a reaction with large heat of activation, 
the velocity of which is $ ~ a“, the flame velocity computed by 
the elementary theory (59) for Dcp / k must be multiplied by 

(Dcp/k) in order to obtain the true value; for , an arbitrary 
function of a must be substituted in $ 

k Hg - H He “ H 

^ ^0 ^ Hg - Bq Hg - Hq ^ ^ 

From the form of the distribution of a and T in space, conclusions 
were drawn by the author concerning the instability of the plane wave 
front relative to its' curving for Dcp^ k (reference 15). liie 
physical properties of such mixtures were confirmed in the experimental 
work of Drozdov and Zeldovich (reference 17). 


8. Equation of An Exothermal Chain Reaction 

The case of several chemical reactions is very common. To de- 
scribe the combustion of a mixture with chain kinetics, considering 
at least two concentrations is necessary, namely those of the initial 
substance a and the active centers and of the active intermediate 
product b. In a number of papers by W. N. Semenov, his followers, 
and the Institute of Chemical Physics a large number of active centers 
were found to be formed in oxidation reactions, thus providing addi- 
tional basis for expecting such to be the case in flames. 

In certain cases, because of definite relations between the con- 
centrations, the problem can be reduced to the case of a single reac- 
tion. Thus, in the case of ozone in the presence of a rapid rever- 
sible reaction, as considered by Le'vrLs and von Elbe (reference 9), 
a balanced thermodynamic relation is established between the concen- 
trations of the active centers and the initial substance. In a 
rapid transfonnation of the initial substance into active centers 
with the liberation of a large amount of heat, the velocity of the 
flame propagation is determined by the single first reaction; the 
succeeding slow transformation of the centers into the final sub- 
stances does not affect the flame velocity. From the mathematical 
point of view, these cases are trivial. In the most frequently 
encountered nontrivial case, however, the transformation of a into 
b generally occurs with a small liberation of heat with a large heat 
of activation and autocatalysis of b; the transformation of b 
into the final products is accompanied by the liberation of a large 
amount of heat and depends little on the temperature. The system of 
equations for an idealized typical chain reaction is used as an 
example . 
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Let the reaction he such that in a closed and heat-insulated 
vessel, the equations of kinetics and the changes in temperatiire 
have the form 


pf =-BiaheE/RT 




dh 

dt 


= Bj^ahe 


-e/rt 


B013' 


> 


dT . ^ ,2 
dt “ ^2^ 


(78) 


The equations of propagation for Dcp » k> y = 0, and for the addi- 
tional simplifying assumptions, k = constant and c = constant, in 
nondimensional variables may be vrritten as 


where 


dz 


dz dz 


dz ^ dz 


B 


+ my - = 0 

(79) 

my ~ - Bpoe”"^!^ «= 0 

(80) 

II - + Bpoe"-'^!^ = 0 

(81) 

. % g-E/RT 
®2 


T 2 - To 

^ rt| 


P 

II 


^■0 
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where aQ is the initial concentration of the substance and the 
remaining nondimens ional variables are the same as previously defined 
with the following conditions: z =» 0, a = 0, 8 = 0, y = 0, ydp/dz = 0, 


z = 1 
a = 1 
P = 0 
y = 0 


do. „ 




The following equation is easily obtained 


d /da. dp t\. /da^dp 
^ dz ^ \ dz dz ” ) ^ \ dz dz ” / '* 


(82) 


where 


a = z - p 

In this manner the two following equations are derived; 


y 


dz 


+ my - p^ 


= 0 


y y + my + BP (z - p) e’-^l^ = O (83) 

dz dz dz 

For 1< z»1/Aj^, set = 0. 

At the present time, the method of finding the value of m in 
such a system is not clear. In the most practical case of the com- 
bustion of hydrogen with oxygen (in this case, all the constants 
entering the equation may be assimied as loiown) , the problem is made 
still more complicated by the fact that Dcp / k and that ^ and 

are introduced in equations (83) so that it is necessary to con- 
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sider a system of three eq,uations.' The principal difficulty, how- 
ever, is in the behavior of the solutions in the neighborhood of 
singular points and the conditions of the simultaneous satisfying 
of the boundary conditions by all the variables . The numerical 
integration of a system of ordinary differential equations for 
different values of one parameter does not offer special difficulties. 


SUMMARY 

1. Equations are set up for heat conductivity and diffusion 
considering the chemical reaction in the flame propagated in an 
explosive mixture. 

2. General properties of the temperatvire distribution in the 
flame are presented. 

3. Conditions are described that must be satisfied by the 
chemical reaction so that a regime of flame propagation with constant 
velocity would be possible. 

4. For a definite relation between diffusion and heat conduc- 
tion, a method is given for determining the flame velocity from the 
equation by numerical integration; the uniqueness of the method 

is also shown. 

5. An estimate in quadratures is given of the upper and lower 
limits of the flame velocity; with Increasing temperatxire dependence 
of the velocity of reaction, both limits approach each other, tending 
toward a value that was given in previous papers on the basis of 
physical considerations. 

6. A method of solution and the limiting value of velocity 
for an arbitrary relation between diffusion and heat conduction 
are given. 

7. Equations are set up for the flame propagation in a chain 
reaction. 
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